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gravitational Chern-Simons term in the action. These are the holographic implementations 
of the usual chiral and the mixed gauge-gravitational anomalies in four dimensional field 
theories with chiral fermions. We discuss the holographic renormalization and show that 
the gauge-gravitational Chern-Simons term does not induce new divergences. In order to 
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^3 , cancel contributions from the extrinsic curvature at a boundary at finite distance a new 

type of counterterm has to be added however. This counterterm can also serve to make the 

Dirichlet problem well defined in case the gauge field strength vanishes on the boundary. 



f^ ' A charged asymptotically AdS black hole is a solution to the theory and as an application 

we compute the chiral magnetic and chiral vortical conductivities via Kubo formulas. We 
find that the characteristic term proportional to T^ is present also at strong coupling and 

- , , , that its numerical value is not renormalized compared to the weak coupling result. 

rS 
H 



5^ I Keywords: Gauge-gravity correspondence, Anomalies, Chern-Simons terms, Transport 



theory 



* karl.landstemer@csic.es 
teugenio.megias@csic.es 
■'■luis.melgar@estudiante.uam.es 
^ fran.penna@uam.es 



Contents 

|l]. Introduction ^ 

§. Holographic Model § 

|3|. Holographic Renormalization 

^ Kubo formulas, anomalies and chiral vortical conductivity ^ 

|5|. Discussion and conclusion 14 



|A] . Codazzi form of Equations of Motion 15 



1^. Technical details on Holographic Renormalization 15 



1^. Equations of motion for the shear sector 20 



|d| . Solutions at zero frequency 21 



1. Introduction 

Anomalies belong to the most interesting and most subtle properties of relativistic quan- 
tum field theories. They are responsible for the breakdown of a classical symmetry due 
to quantum effects. The Adler-Bardeen non-renormalization theorem guarantees that this 
breakdown is saturated at the one-loop level. Therefore the presence of anomalies can 
be determined through simple algebraic criteria on the representations under which the 
chiral fermions of a particular theory transform. In vacuum the anomaly appears as the 
non-conservation of a classically conserved current in a triangle diagram with two addi- 
tional currents. In four dimension two types of anomalies can be distinguished according to 
whether only spin one currents appear in the triangle [1,2] or if also the energy- momentum 
tensor participates [3,4]. We will call the first type of anomalies simply chiral anoma- 
lies and the second type gravitational anomalies. To be precise, in four dimension we 
should actually talk of mixed gauge-gravitational anomalies since triangle diagrams with 
only energy-momentum insertions are perfectly conserved (see e.g. [5]). In a basis of only 
right-handed fermions transforming under a symmetry generated by Ta the presence of 
chiral anomalies is detected by the non-vanishing of dABC = \'~^{Ta{Tb-,Tc}) whereas the 
presence of a gravitational anomaly is detected by the non- vanishing of 6a = Tr(Tyi). 



It recently has been emphasized how at finite temperature and density anomahes 
give rise to new non-dissipative transport phenomena in the hydrodynamics of charged 
relativistic fluids [6-9]. More precisely magnetic fields and vortices in the fluid induce 
currents via the so-called chiral magnetic and chiral vortical conductivities. Although 
there have been many early precursors that found manifestations of this phenomena in the 
physics of neutrinos [10-13], the early universe [14], condensed matter systems [15], the 
recent surge of interest is clearly related to the physics of the quark gluon plasma. It has 
been suggested that the observed charge separation in heavy ion collisions is related to 
a particular manifestation of these anomalous transport phenomena, the chiral magnetic 
effect [16]. The latter describes how a usual (i.e. electro-magnetic) B-field induces via the 
axial anomaly an electric current parallel to the magnetic field. The first application of 
holography to the anomalous hydrodynamics is [17] where the anomalous transport effects 
due to R-charge magnetic fields have been examined. Later studies showed that there is 
also a related vortical effect [7,8], i.e. a vortex in the fluid induces a current parallel to 
the axial vorticity vector w'^ = e^^P^u^dpUx, and related effects of the presence of angular 
momentum had been discussed before in a purely field theoretical setup in [12] and [18]. 
Studies of the chiral magnetic effect using holography have appeared in [19-27] and using 
lattice field theory in [28-30]. A related effect is the so called chiral separation effect that 
induces an axial current in a magnetic field [31]. Chiral magnetic waves have been shown 
to arise through the interplay of chiral magnetic effect and chiral separation effect in [32] . 
Experimental signatures of anomalous transport beyond the charge separation effect have 
been proposed in [33-36]. The experimental status of the observed charge separation in 
heavy ion collision is discussed in [37, 38] . 

A first principals approach to transport theory is via Kubo formulas. In general the 
reaction of a system to external perturbations can be studied via linear response theory. 
The basic objects of linear response theory are the retarded Green functions. Hydrody- 
namic transport coefficients can be extracted from the long-wavelength and low-frequency 
limits of the retarded Green functions. A typical example is the Kubo formula for the shear 
viscosity 

i]=\\m-{TxyT.^y){uj,p = d). (1.1) 

i^->0 id 

Electric or thermal conductivities can be calculated in a similar fashion. 

Also the transport coefficients related to the presence of anomalies can be computed 
via Kubo formulas. For the chiral magnetic effect this has been done in [39] and has been 
further studied in [40]. The Kubo formula for the chiral vortical conductivity was derived 
in [41]. These formulas can easily be generalized to the case of a general non-abelian 
symmetry group, generated by matrices Ta- For the chiral magnetic conductivities and 
chiral vortical conductivities they are 

a% = hm -^ Y. ^-^- (JaJb) (^ = 0,p1 , (1.2) 

a,b 

<y\ = ^lim ^ Yl '-be {jlT'") (^ = 0,p), (1.3) 

a.b 



with a,b,c = x,y, z. 

The properties of these Kubo formulas for the anomalous conductivities are rather 
different from the ones for the dissipative conductivities. Whereas the Kubo formulas for 
the usual transport coefficients are given by the derivatives with respect to the frequency 
at zero momentum, the ones for the anomaly related transport coefficients are given by 
derivatives in momentum space at zero frequency. For this reason the anomalous conduc- 
tivities are contained in the hermitian part of the correlators. In contrast the dissipative 
transport coefficients are contained in the spectral densities, i.e. the anti-hermitian parts 
of the correlators. In the presence of external sources fi{io) that couple to operators O the 
rate of dissipation is describe in terms of the spectral density p^"^ , i.e. the anti-hermitian 
part of the retarded correlator {0^0'^){uj,p), by [42] 

^ = ic./,(-u;)/''H/jH. (1-4) 



This shows that the anomaly related currents due to the conductivities ( |1.2| ) and ( |1.3| ) do 
no work on the system and are therefore examples of dissipationless transport. Let us also 
note that dissipative transport breaks time reversal invariance T whereas anomaly induced 
dissipationless transport preserves T^. 

In [44] these general Kubo formulas were evaluated for a theory of free chiral fermions. 
The results showed a somewhat surprising appearance of the anomaly coefficient bA_ for the 
gravitational anomaly. More precisely the chiral vortical conductivity for the symmetry 
generated by Ta was found to have two contributions, one depending only on the chemical 
potentials and proportional to the axial anomaly coefficient (Iabc and a second one with 
a characteristic T^ temperature dependence proportional to the gravitational anomaly 
coefficient 6^. At weak coupling the anomalous magnetic and vortical conductivities were 
found to be 

B dABC C n r.\ 

^AB = ^^/^ , (1-5) 

V dABC B C , ^Arr2 n an 

^^ = ^^^ ^ + 24^ • ^^-^^ 

This characteristic T^ behavior had appeared already previously in neutrino physics 
[10-13]. It furthermore shows up via undetermined integration constants in effective field 
theory inspired approaches to hydrodynamics. A purely hydrodynamic approach to anoma- 
lous transport has been initiated in [9] (see also [45]). There it was shown that the anoma- 
lous transport coefficient can be fixed by construction of an appropriate entropy current 
whose divergence is positive definite. Later in [46] it was shown that there are additional 
integration constants proportional to T^ and T^. CPT invariance forbids however the T^ 
terms [47] so that at least in systems that can be described by local quantum field theories 
these terms are absent. Recently these studies have been extended to superfiuids [47-49], to 
second order hydrodynamics [43] and to higher dimension [50] where similar undetermined 
integration constants were found. 



""^This point has recently also been emphasized in [43]. 



The usage of Kubo formulas has here a clear advantage, it fixes all integration constants 
automatically. In this way it was possible in [44] to show that the coefficient in front of the 
T^ term in the chiral vortical conductivity is essentially given by the gravitational anomaly 
coefficient b^. The disadvantage of Kubo formulas is of course that we have to calculate 
the potentially complicated correlations functions of a quantum field theory. They are 
easy to evaluate only in certain limits, such as the weak coupling limit considered in [44]. 
In principle the results obtained in this limit can suffer renormalization due to the model 
dependent interactions. 

The gauge-gravity correspondence [51-54] makes also the strong coupling limit easily 
accessible. Strongly coupled non-abelian gauge theories can be described via their gravity 
duals, more precisely in the large A^ and infinite 't-Hooft coupling limit, g'^N — )• oo, allows 
a weakly coupled gravitational description. The drawback is that only some special and 
supersymmetric gauge theories, such as the maximally super symmetric A^ = 4 Yang-Mills 
theory have well understood gravity duals. 

We would like to understand the effects anomalies have on the transport properties 
of relativistic fluids. Anomalies are very robust features of quantum field theories and 
do not depend on the details of the interactions. Therefore a rather general model that 
implements the correct anomaly structure in the gauge-gravity setup is sufficient for our 
purpose even without specifying in detail to which gauge theory it corresponds to. Our 
approach will therefore be a "bottom up" approach in which we simply add appropriate 
Chern-Simons terms that reproduce the relevant anomalies to the Einstein-Maxwell theory 
in five dimensions with negative cosmological constant. ^ 

In this paper we will introduce a model that allows for a holographic implementation 
of the mixed gauge-gravitational anomaly via a mixed gauge-gravitational Chern-Simons 
term of the form 

SCS = I d'x 7^6^^^««^Mi?^ ENPR"" AQR • (1.7) 

Gravity in four dimensions augmented by a similar term with a scalar field instead 
of a vector field has attracted much interest recently [60] (see also the review [61]). A 
four dimensional holographic model with such a term has been shown to give rise to Hall 
viscosity in [62]. The quasinormal modes of this four dimensional model have been studied 
in [63]. 

In section ^ we will define the Lagrangian of our model, derive its equations of motion 
and study how the gravitational anomaly arises. We will find it necessary to add a particular 
boundary counterterm that cancels dependences on the extrinsic curvature of the gauge 
variation of the action. In section ^ we will study the holographic renormalization assuming 
the existence of an asymptotically AdS solution and show that the gravitational Chern- 
Simons terms does not introduce new divergencies. 

In section we will compute the equations of motions for metric and gauge perturba- 
tions in the shear sector. At zero frequency and to lowest order in the momentum these 



^ Very successful holographic bottom up approaches to QCD have been studied recently, either to describe 
non-perturbative phenomenology in the vacuum, see e.g. [55,56], or the strongly coupled plasma [57-59]. 



equations can be solved analytically. This allows us to compute the anomaly related con- 
ductivities in our model. We find that the contribution due to the gravitational anomaly 
is not renormalized compared to the weak coupling result. As is well-known, in a charged 
fluid a current necessarily induces also an energy flux. Kubo formulas for the energy flux 
induced by magnetic flelds and vortices related to energy-momentum correlators. We also 
evaluate these and flnd that no terms proportional to T^ appear. This is of course con- 
sistent with CPT conservation. We also show that the shear-viscosity to entropy ratio is 
unchanged. 

We conclude with a discussion of our results and an outlook towards possible future 
directions in section p. Several technical details of the calculations such as the Gauss- 
Codazzi form of the equations of motion, the details of the holographic renormalization, 
and the equations of motion for the shear sector and their solutions are collected in the 
appendices. 

2. Holographic Model 

In this section we will deflne our model. We start by fixing our conventions. We choose 
the five dimensional metric to be of signature (— , +, +, +, +). The epsilon tensor has to be 
distinguished from the epsilon symbol. The symbol is defined by e{rtxyz) = +1 whereas 
the tensor is defined by e abode = \/—g €{ABCDE). Five dimensional indices are denoted 
with upper case latin letters. We define an outward pointing normal vector n^ oc g -^^ 
to the holographic boundary of an asymptotically AdS space with unit norm hau = 1 so 
that the induced metric takes the form 

hAB = QAB - nATlB ■ (2.1) 

In general a foliation with timelike surfaces defined through r{x) = const can be written 
as 

ds^ = (iV2 + iV^iV^)dr2 + 2NAdx^dr + hABdx^dx^ . (2.2) 

The Christoffel symbols, Riemann tensor and extrinsic curvature are given by 

^NP = ^9^^^ {dNQKP + dpQKM - dKQNp) , (2.3) 

T?A/ _ a -pM Q -pM I -nM -pK -pM -pK (O A\ 

-K NPQ - Opl j^Q- OqI ]^p + l pj^l ]^Q- Iq^l ^p, {ZA) 

Kav = h^Vcnv = -£nhAB , (2.5) 

where £n denotes the Lie derivative in direction of ua ■ 

Finally we can define our model. The action is given by 



S = -^— I d^x./^ 



R + 2K-\fmnF'^''' 



^^MNPQR^^^ (^_FnpFqr + XR^ ENPR"" AQr)\ + SgH + ScSK , (2.6) 

SaH = ^J/.V^hK, (2.7) 

ScsK = -^ j d^xV^XnMe^'''''^''ANKpLDQKJi , (2.8) 



where Sgh is the usual Gibbons-Hawking boundary term and Da = h^V b is the covariant 
derivative on the four dimensional boundary. The second boundary term ScsK is needed 
if we want the model to reproduce the gravitational anomaly at general hypersurface. 
To study the behavior of our model under the relevant gauge and diffeomorphism gauge 
symmetries we note that the action is diffeomorphism invariant. The Chern Simons terms 
are well formed volume forms and as such are diffeomorphism invariant. They do depend 
however explicitly on the gauge connection Am- Under gauge transformations 5 Am = Va/C 
they are therefore invariant only up to a boundary term. We have 

6S = -^ j (i^x^/^^e*^^^^^ [^umFnpFqr + XumR"" enpR"" aqr 



A 

'^^ .la 



[ d''xV^nMe'''^'''^''DMCKpLDQKJi . (2.9) 

Jd 



This is easiest evaluated in Gaussian normal coordinates (see next section) where the metric 
takes the form ds'^ = dr'^ + ^ijdx'^dxK All the terms depending on the extrinsic curvature 
cancel thanks to the contributions from ScsK^- The gauge variation of the action depends 
only on the intrinsic four dimensional curvature of the boundary and is given by 

5S = -^ j d^x^^e'""'^' [^F^nFu + XR' ,„,nR' ^kl) . (2.10) 

This has to be interpreted as the anomalous variation of the effective quantum action of the 
dual field theory. The anomaly is therefore in the form of the consistent anomaly. Since 
we are dealing only with a single U{1) symmetry the (gauge) anomaly is automatically 
expressed in terms of the field strength. We could also express the anomaly in terms of an 
anomalous current conservation equation. One has to be however careful about the defini- 
tion of the current since it is always possible to add a Chern-Simons current and redefine 
J"^ — )• J*" + ce"^"'^''AnFki. This redefined current can not be expressed as the variation of 
a local functional of the fields with respect to the gauge field. In particular the so-called 
covariant form of the anomaly differs precisely in such a redefinition of the current^. A 
good general reference for anomalies is Bertlmann's book [64] where the consistent form of 
the anomaly for chiral fermions transforming under a symmetry group generated by T^ is 
quoted as 



DmJT = mi^e''''T^ 



247r2 



TAdi { AjdkAi + ^AjAkAi 



(2.11) 



with rjH = ± for H € {R, L} for right-handed and left-handed fermions respectively. We 
use this to fix k to the anomaly coefficient for a single chiral fermion transforming under 
a U{1)l symmetry. To do so we simply set Ta = 1 in ( |2.11 ) which fixes the anomaly 
coefficient dABC = ^Tr(T^{TB,Tc}) = 1 and therefore 

' (2.12) 



487rG 967r2 ' 



^Note that the effective field thoery hydrodynamic approaches following [9] typically use the covariant 
form of the anomaly [46] . 



Similarly we can fix A by matching to the gravitational anomaly of a single left-handed 
fermion 

J^ni'J ^ TfiQ 2 ^ nij-ti mkl ; V^-^"^j 

and find 

~167rG ~ 7687r2 ' ^ ' '' 

As a side remark we note that the gravitational anomaly could in principle also be shifted 
into the difFeomorphism sector. This can be done by adding an additional (Bardeen like) 
boundary counterterm to the action 

Sct= fd^xV^A^I^, (2.15) 

with /- = e^'^^'ir^.gdktl+lt^^ptl^tl) fulfilling DrJ"' = \e'^'''R^ nuBJ" rnki- Since this 
term depends explicitly on the four dimensional Christoffel connection it breaks diffeomor- 
phism invariance. 

The bulk equations of motion are 

Gmn - ^9MN = -FmlFn ^ - -F^9MN + 2\eLpQR^M^B {F^^R^ n) '^^) , (2.16) 
V^F^^ = -e^^PQR {kF^pFqu + XR^ enpR"" aqr) , (2.17) 

and they are gauge and diffeomorphism covariant. We note that keeping all boundary 
terms in the variations that lead to the bulk equations of motion we end up with boundary 
terms that contain derivatives of the metric variation normal to the boundary. We will 
discuss this issue in more detail in the next section where we write down the Gauss-Codazzi 
decomposition of the action. 

3. Holographic Renormalization 

In order to go through the steps of the holographic renormalization program within the 
Hamiltonian approach [65,66], first of all we establish some notations. Without loss of 
generality we choose a gauge with vanishing shift vector Nji = 0, lapse N = 1 and Ar = 0. 
So we can use four dimensional (boundary) indices and denote them by small latin letters. 
We therefore also write e{txyz) = +1 and eij^i = V—he{ijkl). In this gauge the bulk 
metric can be written as 

ds"^ = dr^ + jijdx'dx^ . (3.1) 

The non vanishing Christoffel symbols are 

-Hj = Kij = -jij , (3.2) 

rir = Kl (3.3) 

and r*-^ are four dimensional Christoffel symbols computed with jij. Dot denotes differenti- 
ation respect r. All other components of the extrinsic curvature vanish, i.e. Krr = K^i = 0. 



Another useful table of formulas is 

U _ 

J- fcj - 



irj 



-'T' ijk 
-'T' kri 
R^ jkl 



DkK\ + DiKi - 



D'Kki 



-K 



KJ'Ki 



DkKij — DjKik , 



(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 



Note that indices are now raised and lowered with 7jj, e.g. K = Y''^iji s-^d intrinsic 
four dimensional curvature quantities are denoted with a hat, so B} jki is the intrinsic four 
dimensional Riemann tensor on the r = const surface. Finally the Ricci scalar is 



R = R-2K -K^ 



K^jK'^ . 



(3.10) 



Now we can calculate the off shell action. It is useful to divide it up in three terms. The 
first one is the usual gravitational bulk and gauge terms with the usual Gibbons-Hawking 
term. After some computations we get 



5" 



^cs 



n2 



1 



IGttG 

K 



UttG 

8A 
levrG 



d^Xy 



R + 2K + K^ 



K^jK'^ 



:E,E' 






-7e 



ijlk 



AiR^ mklDnK™" + EiKjmDkKj^ + -FikKjmKf 



(3.11) 

(3.12) 

' ;3.13) 



We have used implicitly here the gauge ^^ = and denoted Ai = Ei. The purely four 
dimensional field strength is denoted with a hat. 

Of particular concern is the last term in S^^ which contains explicitly the normal 
derivative of the extrinsic curvature Kij. For this reason the field equations will be gener- 
ically of third order in r-derivatives and that means that we can not define a well-posed 
Dirichlet problem by fixing the 7jj and Kij alone but generically we would need to fix 
also Kij. Having applications to holography in mind we can however impose the boundary 
condition that the metric has an asymptotically AdS expansion of the form 



lij 



,2r 



(0) , - 



■^'^^e- 



4r I 



^+^^) 



H 



(3.14) 



Using the on-shell expansion of Kij obtained in the appendix |B| we can show that the 
last term in the action does not contribute in the limit r — )• oo. Therefore the boundary 
action depends only on the boundary metric 7jj but not on the derivative 7jj. This is 
important because otherwise the dual theory would have additional operators that are 
sourced by the derivative. Similar issues have arisen before in the holographic theory of 
purely gravitational anomalies of two dimensional field theories [62,67,68]. Alternatively 
one could restrict the field space to configurations with vanishing gauge field strength on 
the boundary. Then the last term in Sq^ is absent. We note that the simple form of the 



higher derivative terms arises only if we include ScsK in the action. An analogous term 
in four dimensional Chern-Simons gravity has been considered before in [69]. 

The renormalization procedure follows from an expansion of the four dimensional quan- 
tities in eigenfunctions of the dilatation operator 



Sd 



d xjij 



hv 



(3.15) 



We explain in much details the renormalization in appendix p. The result one gets for the 
counterterm coming from the regularization of the boundary action is 



^ct 



where 



(d-l) 

SttG 


f d^xV-1 
Jd 


h(^-2)^ 








1 / 


pM-p"- l^m ^Ao) '' 


/ J 




4(d-l) iv 


P= - 


R .._ 1 


R) - PS] 




2(d-l)' 


^' (d-2) 





(3.16) 



(3.17) 



As a remarkable fact there is no contribution in the counterterm coming from the gauge- 
gravitational Chern-Simons term. This has also been derived in [70] in a similar model 
that does however not contain ScsK- 



4. Kubo formulas, anomalies and chiral vortical conductivity 

We are now going to evaluate the Kubo formulas for anomalous transport in our holographic 
model. First we note that in a charged fluid a charge current is always accompanied by an 
energy current through 5T^^ = fiSJ"". Therefore charge transport is always accompanied 
by energy transport. Kubo formulas for the energy transport coefficients can easily be 
obtained as well. In [41] it was shown that the chiral vortical conductivity for charge and 
energy transport can be obtained from the retarded Green functions 



|aj=0 



a,b 

a^= lim-^^6,,,(rO«TO^)Uo, 



(4.1) 
(4.2) 



where J* is the (anomalous) current and T*-' is the energy momentum tensor, ay the 
chiral vortical conductivity and ay the vortical conductivity of energy current. The chiral 
magnetic conductivity as and the magnetic conductivity for energy current o"^ are given 
by 



as = lim —-^eabc{J"'J^)\ijj=Q,Ao-- 



a,b 



<'B = i™„^E'-(^°"-'')i 



a;=0 



(4.3) 

(4.4) 



As explained in [23,41] we also have to set the background value of the temporal component 
of the gauge field to zero. Hydrodynamic constitutive relations depend however on a 
particular definition of the fluid velocity. In the case of the anomalous conductivities this 
frame dependence has been addressed in [41] where it was shown how the Landau frame 
conductivities used by Son & Surowka [9] can be obtained from a combination of the 
charge and energy transport coefficient. This combination emerges because of the change 
of coordinates from the laboratory rest frame to a local comoving frame on a element of 
fluid in which there is no energy flux. Applying this change of frame we arrive to the 
transport coefficients in Landau frame 



a,b 



a,b 



(4.5) 
cj=o,ylo=o 



(4.6) 

cj=0 



The frame dependence has also recently been discussed in [50]. The relevant parts of the 
hydrodynamic constitutive relations are 

,5T"" = a% {u"'B'' + u"5™) + a'y (u™w" + u'^uj'^) , (4.7) 

5J" = ubB"^ + ayuj"" , (4.8) 

whereas in Landau frame demanding Um^T"^^^ = we have no contribution to the energy 
momentum tensor but instead 

<5J™ = ^bS" + eyc^™ , (4.9) 

where 5"^ = ie'""'='n„Ffci. 

The AdS/CFT dictionary tells us how to compute the retarded propagators [71,72]. 
Since we are interested in the linear response limit, we split the metric and gauge field into 
a background part and a linear perturbation, 

gMN = OmN + ^ ^MN , (4.10) 

Am = A^^^+eaM. (4.11) 

Inserting these fluctuations-background fields in the action and expanding up to second 
order in e we can read the second order action which is needed to get the desired propa- 
gators [73]. If we construct a vector $^ with the components of qm and Hmn and Fourier 
transforming it 

(27r)« 
it is possible to write the complete second order action on-shell as a boundary term 



^'(^. ^'') = I ^,<ir)e-^'^^'' , (4-12) 



/d A" 
{<^l,AiMJ + ^l,Bjj^i} , (4.13) 

(ivrj" r-5>oo 

where derivatives are taken with respect to the radial coordinate. 
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Now we can compute the holographic response functions from ( 4.13| ) by applying the 



prescription of [71-74]. For a coupled system the holographic computation of the correlators 
consists in finding a maximal set of linearly independent solutions that satisfy infalling 
boundary conditions on the horizon and that source a single operator at the AdS boundary. 
To do so we can construct a matrix of solutions F^ j{k,r) such that each of its columns 
corresponds to one of the independent solutions and normalize it to the unit matrix at the 
boundary. Therefore, given a set of boundary values for the perturbations, 99^, the bulk 
solutions are 

^i{T) = F' j{k,T)^i. (4.14) 

Finally using this decomposition we obtain the matrix of retarded Green functions 

Gij{k) = -2 hm {Aim{F^ j{k, r))' + Bij) . (4.15) 



The system of equations ( 2.16| )- p.l7 ) admit the following exact background AdS 



Reissner-Nordstrom black-brane solution 



d.^ = ^ i-mdt' + dx^) + ^dr-2 , 
^(0) = (p(f)dt = f /3 - ^] dt , (4.16) 



where the horizon of the black hole is located at r = rn and the blackening factor of the 
metric is 

m^i-^.^. (4.17, 

The parameters M and Q of the RN black hole are related to the chemical potential // and 
the horizon fn by 

M = l + | , Q = ^. (4.18) 

The Hawking temperature is given in terms of these black hole parameters as 

^ = Z^/^"H = ^ \ -5 ■ 4.19 

The pressure of the gauge theory is P = xeT^GL^ ^^"^ ^^^ energy density is e = 3P due to 
the underlying conformal symmetry . 

Without loss of generality we consider perturbations of momentum k in the y-direction 
at zero frequency. To study the effect of anomalies we just turned on the shear sector 
(transverse momentum fluctuations) a^ and /i", where a = x,z.^ For convenience we 
redefine new parameters and radial coordinate 

- 4/iAL3 4/iKL3 ^2^2 ^2 



^Since we are in the zero frequency case the fields h" completely decouple of the system and take a 
constant value, see appendix O. 
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Now the horizon sits at u = 1 and the AdS boundary at u = 0. Finally we can write 
the system of differential equations for the shear sector, that consists on four second order 
equations. Since we are interested in computing correlators at hydrodynamics regime, we 
will solve the system up to first order in k. The reduced system can be written as 



= /i? 



u 



hfiu) 



3auB'^{u) + iXkeaji 



(24au3 _ g(i _ ^(^))) 



Bp{u) 



u 



P'i 



+(9aw'^ - 6(1 - f{u)))B'p{u) + 2u(u/if {u))' 
= Bl{u) + ^-i-{B'M) ' ^ ^ 



fin) 



fin) 



+ikeaf 



-^,-x(-ifiu)-l) + 3uAhfiu) + R^ 
uf{u) \a J f{u) 



(4.21) 



(4.22) 



with the gauge field redefined as Ba = aa/lJ-- The complete system of equations depend- 
ing on frequency and momentum is showed in appendix O. This system consists of six 
dynamical equations and two constraints. 

In order to get solutions at first order in momentum we expand the fields in the 
dimensionless momentum p = k/A'KT such as 



B^{u)=Bi^){u)+pB^^\u). 



(4.23) 
(4.24) 



The relevant physical boundary conditions on fields are: /i"(0) = -fT", i?a(0) = B^, where 
the 'tilde' parameters are the sources of the boundary operators. The second condition 
compatible with the ingoing one at the horizon is regularity for the gauge field and vanishing 
for the metric fluctuation [41]. 

After solving the system perturbatively (see appendix ^ for solutions), we can go back 
to the formula ( [4.15| ) and compute the corresponding holographic Green functions. If we 
consider the vector of fields to be 



^l{u) = [B,{u),h^^{u),B,{u),h\{u 
the A and B matrices for that setup take the following form 



(4.25) 



A 



— -^— ^Diag ( -3a/, -, -3a/, - 



(4.26) 



^AdSA 



IQttGL^ 






-AKikp?<j>L^ 
3^5 



-3a 
_ 3 






34 







\ 


-3a 
3 / 



(4.27) 
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BcT 



levrGLS 



/O \ 






3 

















- 


3 

9 r-r 



(4.28) 



u^^/T J 



where B = BAdS+d + ^ScT- Notice that there is no contribution to the matrices coming from 
the Chern-Simons gravity part, the corresponding contributions vanish at the boundary. 
These matrices and the perturbative solutions are the ingredients to compute the matrix 
of propagators. Undoing the vector field redefinition introduced in ( [4.21| ) and ( [4.22| ) the 
non-vanishing retarded correlation functions at zero frequency are then 



Gx,Z 

(^x,tz 
Gtx,tz 



G 



z,tz 



_V3Q_ 

AttGL^ 



i \/^kQ K ik P K 



Gtx,z - 
Gtz,tz 



^z,tx — ^tz,x 

M 



SikQ'^K 2ikXTTT'^ 

+ 



AirGf^ 



G 



167r GL3 



_ iVskQ^K ATxiy/^kQT'^X 

-Gtztx — H Z ^TTH 1- 



27rGf^ 



Gf'n 



(4.29) 
(4.30) 

(4.31) 
(4.32) 
(4.33) 



Using the Kubo formulas ( [4.5^ and 
zero we recover the conductivities 



1) and setting the deformation parameter to 



o-fi 



VSQ K 



27rGf2 



av = <Tb 



"4^2' 



47rG4 



G 



87r2 + 24 



a 



VSQ^K AttVsQT'^X 



M 



V 



27r Gf^ 



Gf^ 



127r2 



+ 



/xr2 

12 



(4.34) 
(4.35) 
(4.36) 



The first expression is in perfect agreement with the literature and the second one shows 
the extra T^ term predicted in [44]. In fact the numerical coefficients coincide precisely 
with the ones obtained in weak coupling. This we take as a strong hint that the anomalous 
conductivities are indeed completely determined by the anomalies and are not renormalized 
beyond one loop. We also point out that the T^ term that appears as undetermined 
integration constant in the hydrodynamic considerations in [49] should make its appearance 
in ay. We do not find any such term which is consistent with the argument that this term 
is absent due to CPT invariance. 
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It is also interesting to write down the vortical and magnetic conductivity as they 
appear in the Landau frame, 



Cb 



iv 




(4.37) 



Finally let us also note that the shear viscosity is not modified by the presence of the 
gravitational anomaly. We know that t] oc limu,_>.o ~ < j'xyrpxy >^^q^ go we should solve 
the system at k = for the fluctuations hi but the anomalous coefficients always appear 



with a momentum k as we can see in ( p. 3D , therefore if we switch off the momentum, the 
system looks precisely as the theory without anomalies. In [75] it has been shown that the 
black hole entropy doesn't depend on the extra mixed Chern-Simons term, therefore the 
shear viscosity entropy ratio remain the same in this model^. 

5. Discussion and conclusion 

We have defined a holographic bottom up model that implements the mixed gauge-gravi- 
tational anomaly in four dimensional field theory via a mixed gauge-gravitational Chern- 
Simons term. We have discussed its holographic renormalization and have shown that the 
Chern-Simons term does not introduce new divergencies. As a first application of this the- 
ory we have computed the anomalous magnetic and vortical conductivities from a charged 
black hole background and have found the T^ terms characteristic for the contribution of 
the gravitational anomaly. 

The most important result is certainly that the numerical values of the conductivities 
coincide precisely with the ones obtained at weak coupling in [44]. This is a strong hint 
towards a non-renormalization theorem for the anomalous conductivities including the 
contributions from the gravitational anomaly. 

We have studied a holographic system with only one anomalous U{1) symmetry. It 
should however present no problem to generalize our calculation to the case with additional 
non-abelian symmetries and various types of mixed anomalies, e.g. mimicking the usual 
interplay of axial and vector symmetries where gauge Bardeen counterterms are necessary 
to implement the correct anomaly structures in the currents [21,23]. 

So far the contributions of the gravitational anomaly have shown up in the calculations 
of Kubo formulas. It is however also possible to calculate directly the constitutive relations 
of the hydrodynamics of anomalous currents via the fluid/gravity approach of [76]. This 
very interesting question can be addressed within the model presented in this paper [77]. 

^For a four dimensional holographic model with gravitational Chern-Simons term and a scalar field this 
has also been shown in [63]. 
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A. Codazzi form of Equations of Motion 



We project the equations of motion ( 2.16| ) and ( |2.17| ) into the boundary surface and the 
orthogonal direction and rewrite them in terms of quantities at the regulated boundary. 
Doing so we get a set of two dynamical equations 

Q = E' + KE' + D,F^' - 4e'^^^ ( kEjFm + iXKWiK.k + 2XR' tkiDsK' 



+4XKksKlDtK' + 4XKstKWiKl , (A.l 



= m + KK] - R) + -E'E. + -F'^^Fj, 



, , 2A + ^E-^Era + -F'™F/, 
(d - 1) V 2 4 



+2A 



2 ' 2 



Ae^'^'-Ef.KfDnK^,^) + le^^^'^D^ (i?fc(i?%)/™ - 2KfK 



^j)m, 



(A.2) 



and three constraints 

-2 



{) = K'^ - Ki-jK'^ -R-2A- -EiE' + -FnF'^ 

2 4 



^t]^ 



+8Ae^^'=' ( Dra{F,jDkKn + Fi.Kkmkr + 2E,KjtDiKi 



k 1 



(A.3) 



= DjK^' - D'K + -EjF^' + 2X€^'"''Dj 



2EkDiKl + Ffcz (^i^^ + KiK'„ 

^^^klmn ^2FklK}DmKi + D, [Fm{R'' nm + 2KKi,) 

+2EkKlR' jni + 2FuKl^{D'Kn^ - DjK'J + 2a,(Ffc/I?ni^i)} , (A.4) 

= D,E' - e'^'^' Upi^Fki + XR' ujR' ski + ^XKi^K^W tki + SXD.KsjDiKl) , (A.5) 
with the notation 



^^v=i(^^•+^/)' 



^ J] •= 2*-^*^' ~^-'' 



(A.6) 



We take Eq. (|A.6| ) as a definition, and it should be applied also when X includes derivatives 
on r, for instance X^^Kij) = \{X^kij + X^kj). 

B. Technical details on Holographic Renormalization 

The renormalization procedure follows from an expansion of the four dimensional quantities 
in eigenfunctions of the dilatation operator 
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5d = 2 I (fx^i~. (B.l) 



This expansion reads 



K] = i^(o) ] + K^2) i + i^(4) i + i^{4) i log e-2^' + • • • , (B.2) 

^i = A(o) i + ^(2) i + i{2) i log e'^*^ + • • • , (B.3) 



where 



(^D-f^(4) j = -4i^(4) ] - 2-ftr(4) j , (50-^(4) J- = -4i^(4) j , 

5d^(0) i = , feA(2) i = —2^(2) i — 2^(2) i , 

6dA(^2) i = -2i(2) i ■ (B.4) 

Given the above expansion of the fields one has to solve the equations of motion in its 
Codazzi form, order by order in a recursive way. To do so one needs to identify the leading 
order in dilatation eigenvalues at which each term contributes. One has 

7ij ~ 0{-2) , Y' ~ 0{2) , E, ~ 0{2) , 4- ~ 0(0) , 

V^~0(-4), Kj~0(0), ff,fci~O(0), Vi~O(0). (B.5) 

Note that for convenience of notation we define 0{n) if the leading eigenvalue of the 
dilatation operator is —n. In practice, in the renormalization procedure one needs to use 
the equations of motion Eqs. (A^) and ( A.3| ) up to 0(2) and 0(4) + 0(4) respectively. Up 
to 0(0) they write 

= Kf,^-K^0))K^0)l-2A, (B.6) 



2A 



= K^)) + K^)K^))--—-5). (B.7) 



Order 0(2) writes 



= 2E:(o)i^(2) - 2K^o) ;-^(2) i-R, (B.^ 

= i^|(2) + ^(0)i^(2) ) + K^2)K^o) )-^. (B.^ 



and finally orders 0(4) and 0(4) for Eq. ( |A.3| ) write respectively 



= 2K(o)K(4) + Kl^ - 2K(o) ;.i^(4) I - K^2) iK(2) i + -F(o) ^Ao) '' ^ (B.IO) 



= 2 (i^(o)i^(4) - i^(o) ;-^{4) 1) log e-^'- . (B.ll) 

The derivative on r can be computed by using 

d 



I d^'xjkm-, = 2 / d^'xE'^jik-^ . (B.12) 

dr J d-fkm J OJkm 
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By inserting in this equation the expansion of -fCj given by Eq. (|B.2| ), one gets d/dr ~ S^ 
at the lowest order. Taking into account this, the computation of Kin\ *■ is trivial if one 
considers the definition of i^jj, i.e. 



K, 



(0) ij 



)7«i|(o) = 2^Dlij = lij 



Then the result up to 0(0) is 



K{o) ] = 5] , 



K, 



(0) 



d. 



(B.13) 



(B.14) 



Inserting this result into Eq. (B.6) or (B.7) one arrives at the well known cosmological 

constant 

d{d-l) 



A 



(B.15) 



We have used in Eq. (|B.7| ) that K^ 
ately from Eqs. (|R|) and (|B.14| ) 



(0) J- = ^dK(q-^ * = 0. The result for K(2) follows immedi- 



K, 



(2) 



P 



R 



2{d-l) 



In order to proceed with the computation of K(2) ) from Eq. 
K^-\/2)- Using the definition of d/dr given by Eq. ( B.12| ), it writes 



(B.16) 
we should evaluate first 



^il(2) 



2 / d'^xK, ' 



(0) milk 



— A-„.+2y...A-„;„-„._A-,„,. 



"• ■^Ikr. 



'5-f, 



-K, 



km 



(2)i 



5dK(2) [ 



-2K, 



mi 



(B.17) 



Because i^(o) ) is the Kronecker's delta, the second term after the first equality is zero, 
while the first one becomes the dilatation operator acting over Kt2) )■ Then one gets from 
Eq. dB^) the result 

Pi = ^J-^ \m - P5il . (B.18) 



K, 



{2)j 



, - , , , R)- P6) 



Note that the trace of K{2) ) agrees with Eq. ( [B.16| ). Using all the results above it is 
straightforward to solve for orders 0(4) and 0(4). From Eqs. ( B.IOJ ) and ( B.ll ) one gets 
respectively 



K, 



1 



(4) 



K, 



(4) 



2{d-l) 
0. 



^^-^'-7^o).,i^o)^^ 



(B.19) 
(B.20) 



In order to compute the counterterm for the on-shell action, besides the equations of 
motion an additional equation is needed. Following Ref. [66] , one can introduce a covariant 
variable 6 and write the on-shell action as 



S. 



1 



on— shell 



SttG 



d*x^/^{K 



(B.21) 
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Then computing Son- shell from Eq. ( B.21| ), and comparing it with the result obtained by 
using Eqs. ( 3.11| )-( |3lT3 ), one gets the following equation 



= e + Ke 



12A 



(d - 1) V 2 ^ 4 ^ J 3 ^ 



(d-1) 



ijkl 



AiR^ mklDnKJ" + EiKjmDkKf^ + -FikKjmK" 



The variable ^ admits also an expansion in eigenfunctions of 5d of the form 

= ^(0) + ^(2) + ^(4) + ^(4) log e-2'^ + • • • , 
where 



(B.22) 



(B.23) 



SD0iO) = 

(4) = - 



2) = 

46'(4) - 261(4) 



-20 



(2) 



fe^(4) 



-40, 



(4)- 



Inserting expansion ( |B.23D into Eq. ( [B.22| ), one gets the following identities 

2A 

= 0(0) +^(0)^(0) --p^r^) ' 

= ^1(2) + -^(2)^(0) + -^(0)^(2) > 

= 0|(4) + i^(4)^(0) + ^(2)^(2) + ^(0)^(4) - TT-l TtF(q) ijF(o) *^ , 



4:{d-iy 



= ^1(4) + ^(0)%) + i^(o)^(4) loge-^^ , 



(B.24) 

(B.25) 
(B.26) 

(B.27) 

(B.28) 



corresponding to orders 0(0), 0(2), 0(4) and 0(4) respectively. Following the same 
procedure as shown in Eqs. ( [B.13D and ( [B.17|) , one gets 



^(0) = ) ^(2) = Sd0(2) = -20(2) 

At this point one can solve Eqs. ( |B.25| ) and ( |B.2(]| ) to get 

P 



'(0) 



^(2) 



(2-d) 



(B.29) 



(B.30) 



Higher orders are a little bit more involved. Using the definition of d/dr, then 0|(4) writes 
0|(4) = 2 / d^xK(o) Llik-^ ^(4) + 2 / d^xK(^4) Llik-^ ^(0) + 2 / d^xK^2) Llik-^ ^(2) 



5d0(4) + 



(2-d) 



d'^xPi 



km 



hkr, 



-P. 



(B.31) 



Note that the second term after the first equality vanishes, while the first one writes in 
terms of 5d- To evaluate the last term at the r.h.s. of eq. ( [B.31 ) we use 



5R = -R'^6jkm + D'D"'6jkm - ^"'DiD'd^, 



km 



(B.32) 
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After a straightforward computation, one gets 



1(4) 



46^(4) - 261(4) + 



{d-l){d-2) 



{d - 2)PjP/ + P' + Di{D'P - D^P]] 



. (B.33) 



Inserting Eq. ( B.33 ) into Eq. ( |B.27| ) one can solve the latter, and the result is 



^(4) 



P}Pi -P'- -/(o) vPm '' + 3 A {D'p - D^p;) 



(B.34) 



The computation of 9\,g^ follows in a similar way, and one gets 9\rg\ = — 46^(4) loge ^^. By 
inserting it into Eq. (B.28), this equation is trivially fulfilled. 

The counterterm of the action can be read out from Eq. ( p. 21 ) by using K and 9 
computed up to order 0(4), i.e. 



S, 



ct 



-s, 



on— shell 



SvrG 



j4 / 

d j;y — 7 



(K(o) - 0(0)) + (i^(2) - 0(2)) + (i^(4) - 0(4)) log e 



-2r 



(B.35) 



From this equation and Eqs. ( |B1^ , ( |B16D , p3oD , ( pISOp and ( p3^ ), one finally gets 



S. 



ct 



jd-l) 
8ttG 



j4 / 

d XV— 7 



1 



1 + 



1 



pip] 

4(d-l) ^ ^ * 



(d-2) 

p2. 



P 



4^(0) *j^(0)'M loge 



-2r 



(B.36) 



The last term in Eq. ( B.34| ) is a total derivative, and so it doesn't contribute to the action. 
As a remarkable fact we find that there is no contribution in the counterterm coming from 
the gauge-gravitational Chern-Simons term. This is because this term only contributes 
at higher orders. Indeed as explained above, in the renormalization procedure we use 
Eqs. (^) and ( p3^ ) up to orders 0(0), 0(2), 0(4) and 0(4), and Eq. ([OD up to orders 
0(0) and 0(2). We have explicitly checked that the A dependence starts contributing at 
0(6) in all these three equations. ^ This means that the gauge-gravitational Chern-Simons 
term does not induce new divergences, and so the renormalization is not modified by it. 



''This result for 6(4) includes a total derivative term which has not been computed in Ref. [66]. To 
compute d{4), in this reference the authors derive the elegant relation ^(4) = ^ ~ ' 7^(4) +7^(4). This identity 
is however valid modulo total derivative terms. 

'^Note that Kj and 9 induce terms proportional to A. Up to order 0(4) + 0(4) these operators write 



^■>\(i)+(i) 



-4K, 



(4) j 



+ 



and e\ 



(4) + (4) 



40(4) + 



where the dots indicate extra terms which are 



A-independent. The only A-dependence could appear in i('(4) * and 6^(4), but these contributions are precisely 
cancelled by other terms in Eqs. ( |A.2| ) and (B.22) respectively, so that these equations become A-dependent 
only at 0(6) and higher. 
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C. Equations of motion for the shear sector 

These are the complete hnearized set of six dynamical equations of motion 



+ik€c 



3 T f 2 , „, , _, ^\ uP' ( - Br{u] 



uf{u) \3a 



Ah^(/(n)-l) + u^ /in«) + '^ 



f{u) 



(C.l) 



= hf {u) - ^^^ - -^^ {k^Kiu) + h^ {u) wk) - SauBM 



U U] 



iXke. 



a/3 



{24au^ - 6(1 - f{u))) ^^ + {9au^ - 6(1 - f{u)))B'i3{u) 



2ulf 



+2u{uhf{u)y - ^ (hl{u)wk + h^iu)k^ 



(C.2) 



= K (u) + 



j/u ^ ujiuY ^ y ' 



uh^;{u) 



+ (9/(u) - 6 + 3an^) 



3\ ""y 



<(^ 



+ 



62 



wkh^{u) + w'^h^{u] 



f{n) ' f{uf 

and two constraints for the fluctuations at ti;, /c 7^ 

= tt; [k^'{u) - 2,auB^{u)\ + f{u)kh^\u) + iAjAe^,^ 2^^ ('ti;/if' + f{u)kh(^'{u) 
+ {9au^-6{l-f{u)))Bi3{u)] . 



(C.3) 



(C.4) 
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D. Solutions at zero frequency 



We write in this appendix the solutions for the system ( [4.2lD -( p:.22 ). These functions 
depend exphcitly on the boundary sources H"^ and Ba, and the anomalous parameters 
R, A. Switching off A we get the same system obtained in [41] 

ikKeap{u — l)a 



h?iu) = H"fiu) 



2(l + 4a)3/2 



(1 + iaf/^u^H^ 



3 Vl + 4au(2au - 1) + 2 (l + m - au^) ArcCoth 



2 + u 



Bk 



kiXeajsiu — 1) 



Bb - 



VI + Aau_ 
ii{u + l)(l + a)TT 3(l + a(5 + a))u {b + 2la + 2a^)- 



2a 



(1 + 4a) 



(1 + 4a) 



3 3i/(w)(l + a(7 + 2a(7 + a)))ArcCoth 
+—i{l + a)TTu — 6au 



2+u 



\/l+4ati 



(u-l)(-l-4a)3/2a 



3/(u)(l + a)^ r 1 21 , 

^ ^^ -—^Log[-l-u + au^] 1 + 



2a(u- 1) 
~gf 2i{u+l){l + afTT 2(l + a)(2 + a(7 + 2a))i 



a2 a(l + 4a) 

(4 + a(25 + a(39 + a(-5 + 4a))))u2 2i{l + afnu'^ 



' + 



a(l +4a) 
u (1 — 5a — 6au) + 

4i/(M)(l + a)(l + 2a)(l + a(5 + a))ArcCoth 



2+u 



^/T+iau 



(u-l)(-l-4a)3/2a2 



2/(u)(l + a)^Log[-l-M + . 
(m - l)a2 

kReap 



Ba{u) = B^+H^u-i 



(h^u (1 + 4a)^/2 + 



2(l + 4a)3/2 
B/3 ( 6aVl + 4au - 2(-2 + a(-2 + 3u))ArcCoth 



2 + u 
Vl + 4au 



+ifcAe 



Q,9 



Bf, 



i(l + a)27r 2(l + a)(l + a(5 + a))u 3i(l + a)7ru 



a(l + 4a) 



-3m 



i(l + a{7 + 2a{7 + a)))(-2 + a(-2 + 3u))ArcCoth 



2a 

2 + M 

\/l+4ati 



(-l-4a)3/2a2 
(1 + a)(-2 + a(-2 + 3u))Log [-1 -u + au^] \ 

2^^2 j 

~p I 4i(l + a)37r (8 + a(48 + a(84 + a(29 + 12a))))u 

^v 3^^ ^ 3a2(l + 4a) ^ 

2i(l + a)^7ru 2(l + a)M2 3 

a^ a 

4i(l + a)(l + 2a)(l + a(5 + a))(-2 + a(-2 + 3u))ArcCoth 



2+u 



^/l+iau 



3(-l-4a)3/2a3 
2(1 + a)2(-2 + a(-2 + 3u))Log [-1 - m + au^] 

3^3 



(D.l) 



(D.2) 
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